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Vibrationally inelastic electron transport 3 ffirdugh "a" m6lecular bridge that is connected to two 



leads is investigated. The study is based on a generic model of vibrational excitation in resonant 
transmission of electrons through a molecular junction. Employing methods from electron-molecule 
. scattering theory, the transmittance through the molecular bridge can be evaluated numerically 

exactly. The current through the junction is obtained approximately using a Landauer-type formula. 
Considering different parameter regimes, which include both the case of a molecular bridge that is 
(N weakly coupled to the leads, resulting in narrow resonance structures, and the opposite case of a 

broad resonance caused by strong interaction with the leads, we investigate the characteristic effects 
of coherent and dissipative vibrational motion on the electron transport. Furthermore, the validity 
of widely used approximations such as the wide-band approximation and the restriction to elastic 
I transport mechanisms is investigated in some detail. 
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I. INTRODUCTION 



The experimental demonstration of the possibility to connect two electrodes by a single molecule and to measure 
a current through such a molecular junction [1-3] has stimulated increasing theoretical efforts to elucidate the basic 
mechanisms of electron transport in such systems (see, for example, [4-8] and references therein). Most of the theo- 
retical work in recent years has been devoted to the determination of the electronic structure of molecular junctions, 
employing a variety of methods that include extended Huckel approaches [9-14], ab initio quantum-chemistry meth- 
ods [15-17], and density functional theory [18-23]. The majority of these studies have focused on elastic mechanisms 
for electron transport, where the current through the molecular junction can be obtained from the single-electron 
transmission probability using the Landauer formula [24-26] . These studies have demonstrated the importance of the 
electronic energy level structure of the molecular bridge: the tunneling of electrons through occupied and unoccupied 
levels of the molecule results in resonance structures in the transmission probability which in turn may cause strongly 
q/Q nonlinear current-voltage characteristics. 

Much less is known about the effect of vibrationally inelastic processes, associated with the vibrational motion of 
the molecular bridge, on the electron transport. In experiments on electron transport through H2 molecules between 
two platinum electrodes [2] as well as Cgo molecules connected to gold electrodes [27], indications for an influence 
of the center-of-mass motion of the respective molecule on the conductivity have been found. Effects of the internal 
vibrational motion of the molecular bridge on the current through the junction, on the other hand, have (to our 
knowledge) not yet been reported. Such effects have, however, been predicted in a variety of theoretical studies. For 
example, the " static" influence of the internal vibrational modes has been studied by averaging the transmittance over 
the probability distribution of the vibrational degrees of freedom [28-30] . The dynamical impact of the vibrational 
degrees of freedom on the tunneling current in molecular junctions has been investigated within, e.g., nearest neighbor 
tight-binding models [31-34]. The effect of the center-of-mass motion on the current has been explored for transversal 
vibrations [35] and for longitudinal vibrations (shuttling mechanism) [36,37]. These studies have demonstrated that 
the vibrational motion of the molecular bridge may result in additional (vibrational) resonance structures in the 
transmission probability which can alter the current-voltage characteristic significantly. Furthermore, the excitation 
of the vibrational degrees of freedom of the molecule provides a mechanism for heating of the molecular junction and 
thus is a possible source of instability [38,39]. 

The effect of vibrationally inelastic processes on electron transport has also been investigated for a variety of closely 
related problems including the tunneling of electrons through long polymer chains (molecular wires) [40,41], electron 
transport through quantum dots and heterostructures [42-46] , as well as the theoretical description of single- molecule 
vibrational spectroscopy in scanning tunneling microscopy (STM) experiments [47-54]. The formally related process 
of electron transport in the presence of a laser field has also been studied [55,56]. 
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Another closely related process is vibrationally inelastic electron-molecule scattering. Here, it is well established 
that the resonant scattering of a low-energy electron from a molecule can result in strong vibrational excitation. This 
process has been studied in great detail experimentally (see, for example, the reviews [57,58]). Furthermore, efficient 
theoretical methods have been developed to describe the interaction of electronic and vibrational degrees of freedom 
in resonant collision processes of low-energy electrons with molecules [59]. As a result, inelastic electron scattering 
from diatomic molecules is now well understood from first principles [60]. Due to the close similarities between the 
process of vibrational excitation induced by electron scattering from a molecule and vibrationally inelastic electron 
transmission trough a molecular junction, it is to be expected that the methods and concepts developed in the former 
field can advantageously be used in the latter field. An example is the so-called wide-band approximation, where the 
energy dependence of the coupling between the molecule and the leads is neglected. This approximation, which has 
been adopted in most of the theoretical studies of vibrationally inelastic electron transport in molecular junctions 
(for exceptions, see Refs. [31,41,43,44,61]) has been tested in detail for resonant electron-molecule scattering [59]. It 
has been shown that the wide-band approximation breaks down at energies close to thresholds and often does not 
accurately describe vibrational excitation processes. 

In this paper we study vibrationally inelastic effects on electron transport through a molecular junction beyond the 
wide-band approximation. To this end, we consider a generic model for vibrational excitation in resonant electron 
transmission processes trough a molecular junction. The model includes the coupling of an electronic resonance state, 
located at the molecular bridge, to the continuum of electronic lead states as well as the coupling of the electronic 
degrees of freedom to a vibrational reaction mode of the molecule. Employing projection-operator techniques [62] 
well-known from electron- molecule scattering [59], the transmission probability through the molecular junction can be 
evaluated numerically exactly within this model. The current trough the bridge is obtained employing a generalized 
Landauer formula [63]. Based on numerical results for models in different parameter regimes, we study the importance 
of inelastic effects on molecular conductance as well as the validity of the wide-band approximation. 

Furthermore, we investigate how vibrationally inelastic effects on electron transport are altered if the vibrational 
motion has dissipative character. Dissipative vibrational processes (such as, for example, vibrational dephasing and 
relaxation) are expected to be of importance in larger molecules or in molecular bridges that are embedded in an 
environment. To describe dissipative vibrational motion, we consider the coupling of the reaction mode to a vibrational 
bath. The various observables are then obtained employing an expansion with respect to the number of quanta in the 
final state of the bath. This technique, which has been proposed some years ago in the context of electron scattering 
from large molecules [64], allows us to describe the effect of vibrational relaxation in an approximate, yet controlled, 
way, without invoking Markov-type approximations. Moreover, it is shown that in the case of identical left and right 
leads and zero bias voltage a unitarity condition can be exploited, which allows a numerically exact evaluation of the 
transmission probability, including the effects of a dissipative vibrational bath. 

This paper is organized as follows: After an introduction of the model and the observables of interest, Sec. II outlines 
the theoretical methods used to describe the transmission probability and the current through the molecular bridge. 
In particular, we discuss various levels of the theoretical treatment: a theoretical description based on purely elastic 
transport mechanisms, the incorporation of vibrationally inelastic processes (coherent and dissipative), as well as the 
wide-band approximation. Sec. Ill presents model studies of vibrationally inelastic electron transport for different 
parameter regimes, comprising both the case of a molecular bridge that is weakly coupled to the leads, resulting in 
narrow resonance structures, and the opposite case of a broad resonance caused by strong interaction with the leads. 
Furthermore, the various levels of theory are critically compared. Finally, Sec. IV gives a summary and concludes. 

II. THEORY 
A. Model Hamiltonian 

To investigate the influence of vibrational motion on the transmission of electrons through a molecular bridge, we 
consider a situation where two metallic leads, which serve as a reservoir of electrons, are connected to a molecule 
through which electrons can be transferred from one lead the other. As has been demonstrated in previous work on 
clastic electron transport (see, for example, [65]), the transmission of electrons through molecular junctions is typically 
characterized by resonances which correspond to the various electronic orbitals of the bridging molecule. From the 
theoretical point of view, the situation is thus characterized by a set of resonance states which are embedded in the 
continuum of lead states. In this paper, we will consider, for simplicity, a situation where only a single electronic 
resonance, corresponding (in the limit of vanishing coupling to the leads) to a molecular anion, contributes to the 
transmission process. 
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From the theory of resonant electron-molecule scattering, it is well known that the influence of vibrational motion on 
the electron transmission can be advantageously described by choosing a basis of diabatic electronic states consisting 
of a discrete state \4>d), which represents the resonance (i.e. the situation where the transmitting electron is situated 
at the bridging molecule) and a set of orthogonal continuum states, \<j>ka), a =L,R, describing the electron in the left 
and right lead, respectively. Accordingly, the Hamiltonian reads 

H S = \<p d )H d ((f> d \ + {\<Pka)(e ka + H )(<p ka \ + \<f>d)Vdka(<f>ka\ + \<Pk a M ka ((f>d\} , (1) 

k,a=L,R 

where H denotes the vibrational Hamiltonian of the neutral molecule in the electronic ground state and Hd the 
vibrational Hamiltonian in the discrete electronic state \(f>d)- The electronic coupling between the leads and the 
molecule is specified by the coupling matrix elements Vdka- 

The electronic parameters of the model Hamiltonian (1) can in principle be determined by electronic structure 
calculations [66]. In the model studies considered below, we have adopted a parameterization which is based on a 
resonance description of the molecular bridge and a simple tight-binding model for the leads, schematically shown in 
Fig. 1. The molecular resonance is described by the discrete state |</> d ) and the states |Z), I = ±1,±2, ... represent the 
atomic sites of the left (— ) and right (+) lead, respectively. It should be noted that in contrast to the lead states, 
the discrete state does not correspond to a single-site tight-binding description (e.g. a single atomic orbital), but 
rather is a molecular resonance state, and thus comprises typically contributions from many atomic orbitals. This 
molecular orbital can be constructed for molecule under bias resulting in voltage-dependent \(f>d}- For information 
on the construction of molecular resonance states we refer the reader to Ref. [59]. The nearest-neighbor coupling 
constants between two lead sites and between the leads and the molecule are specified by [3 and v, respectively, and 
Ml/r denotes the chemical potential in the leads. The stationary continuum states in the right lead are given by 

i« « 2 » 

and similar for the left lead. The energy e of the electron satisfies the dispersion relation 

e = £fca = Ma + 2/3 COS fc, (3) 
in the left and right lead, a — L/R, respectively. Using this particular model for the leads we obtain 



/ sm k 

Vdka = {4>d\H e \\(pka) = Vd ( 4 ) 

As we will see in Sec. II C, the electronic structure of the leads enters the expressions for the observables of interest 
only through the energy-dependent width function of the leads (atomic units with e = H = 1 are used throughout the 
paper unless stated otherwise) 

r Q (e) = 2^]T«5(e-e feQ )|Vd to | 2 . (5) 

k 

The width function r Q (e) is the imaginary part of the self-energy function 

Sa(e) = £ JY^L = Aa(e) _ |r Q (e), (6) 

k 

where e + = e + ry, 7 being a positive infinitesimal. The real part of the self-energy function, the level-shift function 
A a (e), is related to the width function via Hilbert transformation, i.e. 

where P denotes the principal value of the integral. 

For the nearest-neighbor tight-binding model of the leads introduced above, the self-energy function S(z) is given 
by the Hubbard Green's function [67,68], multiplied by the coupling strength v between the last atomic site in the 
leads and the bridge, 
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2v 2 

= =. (8) 

Z - Ma + V ( Z - ^a) 2 - 4 /3 2 

Here, the width of the conduction band is given by 4/3. Analytic continuation in the complex energy (z) plane gives 
the real part 

. , . f^(e-Ma) for |e-Ma| < 2/3 



2/3 2 



(e-M«)T\/(e-Ma) 2 -4/3 2 for ±(e - > 2/3 



and the imaginary part 



r (e) = ! #V4/3 2 -(e-Ma) 2 for |e - Ma | < 2/3 
1 for |e - /i a | > 2/3 



(10) 



As will be demonstrated in Sec. Ill, the inclusion of the energy dependence of the width function r a (e) and thus of 
the coupling matrix elements Vdka, which complicates the theoretical treatment significantly, is crucial in order to 
account correctly for inelastic effects, in particular for energies close to the edge of the conduction band. 

To study vibrationally inelastic effects on the transmission through the molecular bridge, we consider a single 
vibrational (reaction) mode, along which the equilibrium geometry of the discrete electronic state is shifted with 
respect to the continuum states due to the presence of the additional electron at the molecule. Within the harmonic 
approximation, this situation is described by the vibrational Hamiltonians 

A 2 

H =uj s a^a , H d = uj s a^a + \(a + a^) + e d = oj s a d a d + e d . (11) 

Here, u>s is the vibrational frequency of the reaction mode, and and a denote the creation and annihilation operators 
for the reaction mode which are related to the corresponding operators in the discrete electronic state by the shift in 
equilibrium geometry, A/(v / 2ws) , i.e. ad = a + X/los- 

Most of the experimental studies of electron transport through molecular bridges conducted so far have considered 
relatively large molecules with many vibrational degrees of freedom. In large molecules the coupling of reaction 
coordinates (which are strongly coupled to the electronic degrees of freedom) to the remaining (inactive) vibrational 
modes of the molecule results in the process of intramolecular vibrational redistribution, which is well known from the 
spectroscopy of large molecules [69-75]. To study the effect of vibrational relaxation on the electron transmission, we 
adopt a linear response model for vibrational relaxation [64] in the discrete state, where the reaction mode is coupled 
to a bath of harmonic oscillators. Thus the Hamiltonian of the overall system is given by 

H = H S + H B +H SB , (12) 

where the 'system' Hamiltonian H$ is given by Eq. (1), the bath Hamiltonian reads 

H ® =Y^Wjb]bj, (13) 

3 

and the coupling between the reaction coordinate and the bath is given by the bilinear interaction 

H SB = \M ]T Cj (a d b] + afo-NM (14) 

3 

Here, bj and bj denote creation and annihilation operators for the bath mode with frequency uij and Cj is the 
corresponding system-bath coupling constant. There is no bilinear coupling between the system and the bath modes 
if the molecule is in the neutral state (and thus the electron in the lead states) in our form of Hsb- In reality there 
is, of course, such coupling and also relaxation due to anharmonic effects in the neutral state of the molecular bridge, 
describing, for example vibrational relaxation processes after an electron has been transmitted through the bridge. 
The study of these effects, which for the process of resonant electron-molecule scattering were found small compared 
to the relaxation mechanism in the discrete molecular state [64], will be the subject of future work. 

As is well known [76,77], all properties of the vibrational bath which influence the dynamics of the system are 
characterized by the bath spectral density 
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JH=^^( W - Wi ). (15) 

3 

In the numerical calculations reported below, it will be modeled by a continuous Ohmic bath with exponential cutoff 
[76,77] 

J(w)=Tjwe- w / We . (16) 

Here, the characteristic frequency uj c defines the maximum of the spectral density and the overall strength of the 
system-bath coupling is measured by i]. 

B. Observables of Interest 

Various observables are of interest when investigating the influence of vibrational motion on the electron transport 
trough a molecular bridge. The most detailed information on the transmission process of a single electron is comprised 
in the initial- and final-state resolved scattering probability. Employing scattering theory, it is straightforward to show 
that the probability for scattering of an electron with the energy ej from lead on into a state with energy e/ in the 
lead a/, accompanied by a vibrational transition from state |u»)|vj) to state l^/^v/) 1 (we denote the bath vibrations 
with bold face letters and the system mode vibrations with italic letters) is given by [59] 

taf<-a t (£f,Vf,Vf,ei,Vi,Vi) = 5(ej + E w% + E Vi - e f - E Vf - E Vf )T ai (ei)T af (ef) 

x|(v / |( U/ |^ d |(e+-^)-Vd)K)|v i )| 2 ) (17) 

where E Vi , E Vi and E Vf , E Vf are the initial and final vibrational energies of the bath and system modes, respectively. 
Writing the transition probability in the form (17), we have assumed, for simplicity, that the coupling element Vdka 
does not depend on the vibrational degrees of freedom. The generalization of the formula to include such effects is 
straightforward (see, for example, Ref. [59]). 

While Eq. (17) describes the most detailed information on the scattering process, in experiments typically more 
averaged observables are measured. It is thus expedient to introduce the integral transmission probability from lead 
a, into lead a/ (summed over all possible final vibrational states) 

t af ^ ai (e f ,ei) = ^ t af <- ai (ef,Vf,Vf,ei,Vi,\i), (18) 

Vf.Vf 

and the total transmission probability, integrated over the final energy of the electron, e/, 

t af ^ ai (ei) = J t a/ «_ a4 (e/,ei)de/. (19) 

The most important observable for the study of electron transport trough a molecular bridge is, of course, the 
current which is induced when a finite voltage is applied to the molecular junction. To calculate the current through 
the bridge, we employ the generalized Landauer formula [5] 



I= \J de * / de/{tfl^(e / ,e i )/ i (e i )[l-/ fl (e / )]-t i ^ fl (e / ,e i )/ fl (e i )[l-/ i (e / )]} 



(20) 



where f a (E), a = L,R, denotes the Fermi-Dirac distribution for the left and right lead, respectively. 

In principle the basis states \4>d), \<l>ka), and therefore also the functions T a (e) and vibrational Hamiltonian H dl 
depend on the bias voltage V across the bridge. Here, we assume for simplicity that the bias voltage V enters Eq. 
(20) only through the Fermi distribution of the leads and the width function T a (e) [cf. Eq. (3), (5)] via the chemical 
potentials of the leads Hl/r = ±V/2. Moreover, in this paper we will not consider thermal effects on the electron 
transport, i.e. in all numerical calculations reported below we have taken T = K. Thus, the initial state of the 
vibrational degrees of freedom is the ground state of the system and bath modes, respectively, |u, = 0)|vj — Ob) and 
the Fermi-Dirac distribution f a (E) in Eq. (20) reduces to the Heaviside step function. 



1 The formula has to be modified for elastic backscattering, i. e. a; = a/, |ui)|v») = |«/)|v/), confront (31b) 
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The validity of the formula (20) for calculating currents including inelastic effects has been discussed controversially 
in the literature. The current formula is sometimes used without the Pauli exclusion principle factors [1 — /(e)] (see 
for example the discussion in Refs. [26,78,79]). Employing nonequilibrium Greens-function theory and the second 
quantized version of the Hamiltonian (1), it can be shown that Eq. (20) gives the correct description of the current 
in the limit where many-electron processes are negligible for the dynamics [80]. This limit is sometimes called the 
single particle approximation. For the study of the tunnelling through a vibrating single-molecule junction beyond 
this approximation (but within the wide-band limit) see [81] and references therein. It should also be mentioned that 
the use of Eq. (20) implicitly assumes that the bridge relaxes into the ground vibrational state \vi — 0}|vj = 0_b) 
before each subsequent transmission event. In view of the low currents observed experimentally, this should be a 
reasonable assumption (for a similar discussion in the case of STM currents, see, for example, Ref. [50]). 



C. Method of Solution 



We explain the method of solution of the problem in this paragraph, starting from simplified cases which are later 
also discussed and shown for comparison with the final solution. 



1. Elastic transmission 



Let us first consider the case without coupling to the vibrational degrees of freedom [i.e. A = rj = in Eqs. (11), (16)]. 
In this case, only elastic processes contribute to the electron transport, and the total transmission probability is given 
by 

+ ( c \-+ ( c \- r K (e t )r L (e t ) 

R ^ L[ i) L ^ R{ l) h - e d - A L ( £l ) - A R (ei)] 2 + [T^) + r K (e,)] 2 /4 ' ( ' 

which is the well-known result for elastic resonant tunneling transmission. The Landaucr formula (20) for the current 
is exact in this case (see, for example, Ref. [82]) and the evaluation of the current reduces to a simple numerical 
integration. 



2. Inelastic transmission without dissipation 



Let us next consider the influence of the coupling to the vibrational degrees of freedom on the transmission of 
the electron, i.e. vibrationally inelastic transmission through the molecular bridge. If we exclude dissipative processes 
induced by the coupling to the bath, the dynamics is described by the system Hamiltonian Hs, Eq. (1). The probability 
for transmitting an electron with energy q from the left lead into a state with energy e/ in the right lead, accompanied 
by a vibrational transition from the state \vi = 0) (which is the initial vibrational state at temperature T = 0) to the 
state \vf) is given by ( [79,83]) 



tR^Lfa,*) = ^<5(e 4 - e f - E v/ )r R (e / )r L (e i ) (v f \G { d S) (ci)\0) 



v f 



where 



G s d {E) = (4> d \(E+ - Hs)- 1 ^) 



(22) 



(23) 



denotes the Green's function projected on the resonance state. The exact one-particle transmission probability (22) 
can be obtained from the solution of both time-dependent [51,61] and time-independent Schrodingcr equation (see for 
example [43,44,79,83]). The Green's function (23) can also be written in closed form. Employing projection operator 
techniques well known from the theory of resonant electron-molecule scattering (see [59] and references therein), 
G%(E) can be recast in the form 



G s d (E) = \E + -H d -S L {E- H ) - Z R (E - H ) 



(24) 



This form has the advantage that the electronic continuum has been formally eliminated and only vibrational dynamics 
in the discrete electronic space has to be evaluated. We would like to stress that we do not assume the wide-band 



6 



approximation (see the following section) in Eq. (24) and the Hamiltonian operator Ho enters the energy dependence 
of E. Efficient techniques have been developed [84,86] to evaluate the matrix elements of the Green's function (24). 
In the present case, where H and Hd describe harmonic oscillators, the Green's function (24) can be obtained, e.g., 
by inverting a basis representation of the operator E + — Hd — Y,r,(E — H$) — Y,r(E — Ho) for each energy E, employing 
efficient algorithms for the inversion of tridiagonal matrices [87]. 



3. Wide-band approximation 



While the expressions discussed so far take full account of the energy dependence of the width function T(e), in the 
majority of previous work on the effect of vibrational motion on electron transmission, the so-called wide-band (WB) 
approximation has been invoked, where the width function is assumed to be constant, i.e. r(e) = const (for recent 
examples see [85,35]). In this approximation, the level-shift function vanishes [cf. Eq. (7)], A(e) = 0. Introducing the 
eigenstates \n) and eigenenergies E n for the operator Hd the transmission probability, Eq. (22), can be written as 

2 



Cl^ = T R r L ]T 5(e t e f - E Vf ) 



E 



(v f \n)(n\0) 



ti — E n 



(r R + r L ) 



which yields after integration over the final energy of the electron the total transmission probability 

I HO) I 2 



4r L ) (e i ) = r R r L ^ 



r R r L ]T 



(a - E n y + \ (r R + r L )2 



(25) 



(26) 



(a - E n y + \ (r R + r L )2 



In Sec. Ill, we will study the validity of the wide-band approximation based on the comparison of Eq. (26) with the 
full inelastic transmission probability given by Eq. (22). 



4- Inelastic transmission including vibrational relaxation 

Finally, we consider the case of inelastic transmission of an electron through the molecular bridge in the presence 
of vibrational relaxation. The theoretical treatment of this problem is considerably more complicated than the 
cases considered above, because we have to deal with two qualitatively different continua: the electronic scattering 
continuum describing the leads and the dissipative vibrational mode continuum of the bath. In principle, it is possible 
to start from the general formula for the transition probability (17) and derive formal expressions as in Eqs. (22), (24) 
with the bath modes included in Hg and Hd- However, with increasing number of bath modes it becomes difficult to 
perform the matrix inversion in Eq. (24) and, therefore, such an approach is limited to very few bath modes. 

To circumvent this problem, we adopt an approach that has been proposed to describe the effect of vibrational 
relaxation in the context of resonant electron scattering from large molecules [64]. The basic idea of this method is 
to express the transmission probability, Eq. (18), as a sum 

oo 

*R^L(e/, €i) = 4^l( £ / j ( 27 ) 

where t R ^ L (e/,e i ) describes transmission processes with m quanta of excitation in the final state of the bath (as 
was mentioned above, we assume that the initial state of the bath at temperature T = is given by |vi = 0^)). 
In the parameter regime where the vibrational relaxation rate is small compared to the electronic decay rate of the 
resonance, only the first few terms in the expansion will contribute to the overall transmission probability. The first 
two terms in the expansion, which correspond to processes where the final state of the bath contains zero or one 
quantum of excitation, respectively, read 

ilL(^^0=E^- e /-^/) r R( e /) r L^)K w /l G d^)|0)| 2 , (28a) 

v f 

fiu-L - E J ( £ < - e f -K / )r R (e / )r L (e i ) \{v f \G d {e f +E Vj )a d G d {e t ) |0)| 2 , (28b) 

Vf 
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with the Green's function 



G d (E) = (cf, d \ (e+ -H s - \<f> d )a d (<f> d \ J du E+ ^H a _J <l>d)ad(4>d\) I'M (29) 

= {m(e + -H s -\4> d )a\ J dw J(u) G s d (E - w) a d (cj> d \^ \<t>d)- 

The expressions for higher order terms (to > 1) can be found in Ref. [64]. The Green's function G d {E) can be 
evaluated either as the unique solution of the Lippmann-Schwinger equation 



G d (E) = G s d (E)+G s d (E) 



4 J dujj(uj) G s d {E-w)a d 



G d (E) (30) 



using the iterative Schwingcr-Lanczos method [86,88] (which converges with only few iterations) or employing matrix 
inversion techniques as described in Ref. [64] . 

In general situations where the left and right lead are not identical, e.g., due to a non-zero bias voltage, the expansion 
in Eq. (27) has to be terminated for practical reasons at a certain order m and is thus only applicable in the case of 
a weakly damped system mode. We would like to emphasize, though, that the electronic coupling of the molecule to 
the leads is treated exactly in the approach outlined above and it can thus be arbitrarily strong. Furthermore, the 
treatment of the dissipation described by Eq. (28) does not invoke Markov-type approximations and is, therefore, not 
limited to situations where the bath-correlation time is short compared to the system dynamics. 

In the case of identical left and right leads and zero bias, however, one can exploit the unitarity condition to express 
the sum of all higher order corrections, i.e. t^_ L with m > 0, in terms of the elastic (with respect to the bath) terms 

t^-L an d ^lI-l- T° see tms ' ^ i s n °ted that in the symmetric case we have t^_ L — t^_ L for all m > 0. For m=0, 
the total transmission probabilities are given by 

oo 

4°Il(^) = £ r R ( £j - E Vf )Y^) \{v f \G d {em\ 2 , (31a) 

v f =0 

OO 

4°2-Lte) = I 1 " *r L (6,)(0|G d ( £j )|0)| 2 + r L ( £l - E Vf )Y^e t ) |( W/ |G d ( £l )|0)| 2 . (31b) 

Vf = l 

The difference between the transmissions probabilities from the left to the right lead and from the left back to the left 
lead is due to the special role of the transition amplitude with the same initial and final state which only contributes 
in the latter case. Together with the unitarity condition 

*£LlM + E + 4°l L (e) + E 4 m 2i» - 1, (32) 

m>0 m>0 

this yields 

tfl-L(*) = \ (l + t£l L (ei) - 41l(^)) = -r L (eO Im<0|Gd(ci)|0>. (33) 

Eq. (33) is an exact formula for the total transmission of the electron through the molecular bridge, including the bath 
to all orders in the system-bath coupling. Moreover, it can easily be evaluated numerically. Although it is limited 
to symmetric leads, and thus cannot be applied directly to calculate the current, it is very helpful for checking the 
convergence properties of the expansion (27) for zero bias or in the linear response regime. 



III. RESULTS AND DISCUSSION 



In this section we present the results of a model study of the influence of vibrational motion and vibrational 
dissipation on the transmission probability and the current-voltage characteristics of a molecular bridge. To obtain 
a comprehensive picture of the various mechanisms, we shall consider models in different parameter regimes. In 
particular, we will consider both the case of a molecular bridge that is weakly coupled to the leads, resulting in 
narrow resonance structures, and the opposite case of a broad resonance, caused by strong interaction with the leads. 
The parameters specific to the different models are collected in Table I. The nearest-neighbor coupling strength in 
the leads is (3 — 1 eV in all models considered (this parameter is an overall energy scaling factor). The characteristic 
frequency of the bath is chosen to coincide with the frequency of the system mode, i.e. u>c — ws, and we will consider 
a relatively weak coupling between the system mode and the bath, — 0.1. 
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A. Tunneling through a narrow resonance 

First, we consider the resonant transmission of electrons in the tunneling regime, which is characterized by a 
relatively weak coupling of the electronic state localized at the bridge to those in the leads. We have chosen a 
coupling strength of v = 0.2 eV, which corresponds to a fifth of the nearest-neighbor hopping amplitude f3 in the 
leads. We will, furthermore, consider a localized state that is situated well inside the conduction band with an energy 
of 64 = 0.5 eV (for zero voltage, the conduction band extends over the range [ -2 eV, 2 eV]). 

We start with Model A, which is characterized by a relatively weak coupling between the electronic degrees of 
freedom and the vibrational motion of the system mode, A = 0.3 eV. Fig. 2 depicts the transmission probability for 
model A obtained for zero voltage using the different levels of theory introduced in Sec. II C: the elastic transmission 
probability, Eq. (21), the vibrationally inelastic transmission probability integrated over the final electron energy, 
Eq. (22), the wide-band approximation of the vibrationally inelastic transmission probability, Eq. (26), and the 
transmission probability in the presence of vibrational relaxation, Eq. (33). The elastic transmission probability 
(dotted line) exhibits a rather narrow peak at the position of the discrete electronic state. Including the coupling of 
the electronic degrees of freedom to the vibrational motion of system mode (vibrationally inelastic transmission, thick 
dashed line), this peak is seen to become split into several sub-peaks which correspond to the different vibrational levels 
in the discrete electronic state. This effect is well known from previous studies [45,46]. The comparison between the 
clastic and inelastic transmission probability in Fig. 2 demonstrates that even in the case of relatively weak electronic- 
vibrational coupling a theoretical treatment which only includes elastic processes provides a rather poor description. 
The wide-band approximation (thin dashed line), on the other hand, which takes the vibrational excitation into 
account but neglects the energy dependence of the width function, is seen to give an excellent description of the 
transmission probability. This is due to the fact that in model A the resonance is situated well inside the conduction 
band and, therefore, threshold effects are are negligible. 

The inclusion of the coupling to the vibrational bath (solid line) causes a further broadening and slight shift of the 
peaks belonging to the first and second excited vibrational level in the discrete electronic state. The main peak, which 
corresponds to tunneling trough the ground vibrational state of H c i 1 on the other hand, remains almost unaffected. 
This is due to the fact that (in a zero-order picture without electronic coupling) the ground vibrational state of Hd 
is a stationary state, while all excited vibrational states decay into the ground state. 

The result for the inelastic electron transmission including vibrational relaxation depicted in Fig. 2 has been obtained 
employing Eq. (33) which is only valid for zero bias. For situations with non-zero bias (in particular, to evaluate the 
current through the bridge), we will use the expansion of the total transmission probability in terms of the number of 
excitations in the final states of the bath, Eq. (27). It is therefore important to study the validity of this expansion. 
The results in Fig. 3 demonstrate that for the present example the expansion (27) is well converged if bath states 
with zero and one quantum of excitation (to = 0, 1) are taken into account. 

Fig. 4 shows the current through the bridge as a function of the applied voltage for model A. The results have been 
obtained employing Eq. (20). As has been discussed in detail by other authors [89,90], each resonance peak in the 
transmission probability has its counterpart in a step in the current-voltage curve. Thereby, the steps occur in the 
order as they appear in the transmission probability counted from the zero- voltage Fermi energy (which is set to zero). 
Accordingly, the current based on the elastic treatment of the transmission (dotted line) exhibits only a single broad 
maximum corresponding to the position of the discrete electronic state. The inelastic current (thick dashed line), on 
the other hand, exhibits several steps which belong to the different vibrational peaks in the transmission probability. 
The wide-band approximation (thin dashed line) is seen to give very good results except at high voltages, where the 
resonances are closer to the edge of the conduction band and, therefore, the energy dependence of the width function 
becomes important. 

Although there is a pronounced effect of the vibrational motion on the current in model A, the influence of the 
coupling to the bath (thick solid line) is rather small. This is a consequence of the location of the discrete electronic 
state, which is situated 0.5 eV above the Fermi energy in this model. As a result, the current for low voltages (0 — IV) 
is almost exclusively due to tunneling of electrons via the resonance corresponding to the ground state of Hd, which 
is hardly affected by dissipation (cf. the discussion above). 

If we change the discrete state energy to = —0.5 eV (Model B), the transmission functions, depicted in Fig. 5, 
remain virtually unchanged except for a shift in energy by 1 eV. In contrast, the current-voltage characteristic for 
Model B, shown in Fig. 6, is qualitatively different from that of Model A (cf. Fig. 4). The reason is that the order of 
the peaks as counted from the Fermi energy is reversed and therefore the low- voltage region of the current is influenced 
by tunneling trough excited vibrational states of Hd which are more strongly affected by the presence of the bath. 
The current with and without vibrational relaxation thus differs by more then 50 % for voltages in the range 0.5 — 1 
V. 

The differences between the various levels of the theoretical treatment become more significant if wc consider a 
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model with larger coupling between the electronic and vibrational degrees of freedom. Such a situation is described 
by model C, where the vibrational frequency is chosen as luq = 0.4 eV and the vibronic coupling strength as A = 0.7 
cV. The location of the discrete state is the same as in model A, td — 0.5 eV. The transmission probability for model 
C is shown in Fig. 7. It is seen that the stronger vibronic coupling results in a pronounced vibrational progression 
in the transmission probability. In contrast to the cases considered above, the wide-band approximation essentially 
fails to describe the transmission probability. It predicts an amplitude which is too small by about a factor of two. 
Furthermore, the position of the peaks is not correctly described in the wide-band approximation due to the neglect of 
the level-shift function A(E), and the individual peaks in the numerically exact results are narrower than in the wide- 
band approximation. This effect of 'vibrational narrowing' is well-known from resonant electron-molecule scattering 
[87] and results from the interference between overlapping resonances. 

The effect of vibrational relaxation on the transmission probability, illustrated in Fig. 8, also is much more pro- 
nounced in this model than in model A. Except for the lowest two peaks, the vibrational resonances are smeared 
into a broad hump, when vibrational dissipation is included. This is due to the fact that the vibrational relaxation 
process becomes more effective for higher excited vibrational states. Also shown in Fig. 8 is the expansion of the total 
transmission probability in terms of the number of excitations in the final state of the bath, Eq. (27) . As a result of 
the importance of higher vibrational states and the relatively small electronic coupling to the leads, the expansion is 
seen to converge much slower than in the models considered above. 

Let us next consider the current-voltage characteristic for model C, depicted in Fig. 9. It is seen that the coupling to 
the vibrational motion has a rather strong effect on the current through the bridge. In particular, the current-voltage 
characteristic exhibits step-like structures corresponding to the different vibrational levels in the discrete electronic 
state. As expected from the discussion of the transmission probability above, both the elastic treatment and the 
wide-band approximation fail severely in the description of the current. Fig. 10 demonstrates that the coupling to 
the bath has a strong effect on the current-voltage characteristic in model C. In particular, the step-like structure 
is washed-out and the magnitude of the current decreases by more than a factor of two. Also shown in Fig. 10 is 
the contribution of the different terms in the expansion (27) to the total current. In contrast to the transmission 
probability, the expansion for the current converges relatively fast for this model. The reason for this at first sight 
surprising finding is that, due to the Pauli principle exclusion factors in the formula for the current, terms with a 
higher number of excitations in the final state of the bath are suppressed at lower voltage. 

B. Transmission through a bridge strongly coupled to the leads 

All models considered so far were characterized by a relatively weak coupling of the discrete electronic state to the 
leads, which results in narrow resonance structures. In this section we shall consider the opposite limit of a localized 
state that is strongly coupled to the leads. Such a situation is realized in model D, where the coupling strength 
between the discrete electronic state and the leads is chosen as v — 1 eV. The position of the discrete state is in the 
upper part of the conduction band, td = 1-6 cV. All other parameters are the same as in Model C. We mention that 
if the coupling to the bath is not considered, model D is essentially equivalent to a model studied by Gelfand et al. 
[61] in the context of inelastic tunneling in heterostructures. 

The results for the transmission probability are depicted in Fig. 11. The strong coupling to the leads results in a 
rather broad transmission probability, which is qualitatively well described taking into account only elastic processes. 
Vibrationally inelastic contributions to the transmission probability manifest themselves in various cusp structures. 
Thereby, each cusp indicates the opening of a new vibrational channel. It is well-known from the theory of electron- 
molecule scattering that the wide-band approximation is not at all applicable in this case. The effect of the coupling 
to the bath, which is well described including the two lowest term in the expansion (27), is very small. This is a 
consequence of the strong electronic coupling which results in a very short residence time of the electron on the bridge. 

Fig. 12 displays the current-voltage characteristic for model D. As to be expected from the transmission probabilities, 
the differences among the various levels of theory are small with the exception of higher bias voltages. For higher 
voltages, the elastic current vanishes due to the empty overlap of the left and right conduction bands. Inelastic 
transmission processes, however, which are accompanied by an energy loss of the electron, are still allowed. 

Although there are no structures in the transmission probability in Fig. 11 which are obviously related to the 
position of the discrete electronic state, the location of this state does play an important role. This is demonstrated 
in Fig. 13, which shows the transmission probability for Model E, which differs from Model D only in a lower energy 
of the discrete electronic state, = —0.7 eV. This different location of the discrete state results in pronounced peaks 
and minima in the transmission function, which are somewhat smoothed, but not destroyed, by the coupling to the 
bath. To facilitate the interpretation of these structures, Fig. 14 shows the potential-energy curve of the discrete state 
(corresponding to Hj) together with the energies of the vibrational states. In addition, the potential energy of H is 
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shown, shifted by ±2 eV, respectively, to indicate the energy which electrons coming from the conduction band may 
carry into the bridging molecule. Though the localized state is strongly coupled to the continuum, due to the shift of 
the two potential curves, the ground vibrational state in Hd has only a small overlap with the respective ground state 
in H - Consequently the coupling between the two states is effectively small and sharp resonances may be observed. 
If the discrete state is localized higher in energy (as in model D), the potential energy curve of Hd is shifted up. Then, 
the ground state of Hd still has a small overlap with the ground state of H but the decay into higher vibrational 
states becomes energetically possible. The sharp resonances thus "dissolve" in the continuum. As is demonstrated in 
Fig. 15, the sharp structures close to the bottom of the conduction band have no significant effect on the current, and 
the overall appearance of the current- voltage characteristic is similar as in model D. 

IV. CONCLUSIONS 

In this paper we have studied vibrationally inelastic effects on electron transport through a molecular bridge that 
is connected to two metal leads. The study was based on a generic model for vibrational excitation in resonant 
electron transmission processes through a molecular junction. Employing projection-operator methods well-known 
from resonant electron-molecule scattering, we have outlined how the transmission probability can be evaluated 
numerically exactly within this model, without invoking the wide-band approximation or perturbation theory with 
respect to the coupling between the bridging molecule and the leads. Furthermore, the influence of dissipative 
vibrational processes was investigated by considering the coupling of a vibrational reaction mode to a dissipative 
bath. 

The results of the model study can be summarized as follows: In the case of tunneling through a molecular 
bridge which is weakly coupled to the leads, the transfer of an electron may result in strong vibrational excitation, 
which manifests itself in pronounced vibrational resonance structures in the transmission probability and in a step- 
like appearance of the current-voltage characteristic. Since in this case the residence time of the electron on the 
molecular bridge is relatively long, dissipative processes such as vibrational relaxation can have a significant effect 
on the dynamics. In particular, they result in a broadening of the resonance peaks in the transmittance and of the 
step-like structures in the current-voltage characteristic. Furthermore, vibrational relaxation may result in this case 
in a significant reduction of the overall magnitude of the current. Due to the pronounced effects of the vibrational 
degrees of freedom, a theoretical treatment which only includes elastic processes is not appropriate in this parameter 
regime. Our studies also show that the wide-band approximation can only be applied if the electronic resonance state 
is situated well within the conduction band and the electron-vibrational coupling is weak. 

In the opposite case of a molecular bridge that is strongly coupled to the leads, the transmission probability 
is typically characterized by a broad distribution, which in turn results in a rather structureless current-voltage 
characteristic. Nevertheless, the vibrational motion may manifest itself in cusp structures in the transmittance. 
Furthermore, sharp resonance structures may occur in the transmission probability, if the energy of the discrete 
electronic state is low enough and the electronic-vibrational coupling sufficiently strong, such that the vibrational 
ground state of Hd has some overlap with lower- lying vibrational states of H . Except for the latter case, the effect 
of vibrational relaxation is very small in this parameter regime. The comparison of the results obtained at different 
levels of theory shows that methods which only include elastic processes can give a rather good qualitative description 
of the electron transport in this case, although they miss the detailed cusp and resonance structures. The wide-band 
approximation, on the other hand, is not valid in this parameter regime; due to the strong molecule-lead coupling, 
threshold effects become important which are neglected in the wide-band approximation. 

To study the basic mechanisms of vibrationally inelastic electron transport, we have focused in this work on 
relatively simple models with a single harmonic reaction coordinate and a single electronic resonance state. It should 
be noted, however, that the methods employed in this work are not limited to these models. The extension of the 
theory to an anharmonic reaction coordinate, several reaction coordinates, and several resonance states is relatively 
straightforward. Also, it should be emphasized that the potential-energy surfaces of such models can in principle be 
determined ab intio by electronic structure calculations. In this way, for example, the possibility of dissociation of 
the molecular bridge induced by a strong current can be studied. 

Finally, it is noted that in the present work the current through the molecular junction was obtained with the 
generalized Landaucr formula, Eq. (20). Although this formula gives the correct description in the limit of weak 
coupling between molecule and leads as well as in the situation when only single-electron processes are important, it 
needs to be extended for applications where these assumptions are not fulfilled. A theoretical treatment of inelastic 
processes without these limitations is possible within the framework of non-equilibrium Green's function theory [82] . 
The combination of this formalism with the methods employed in this paper is a challenging subject for future research. 
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Model 




V 


u>s 


A 


A 


0.5 


0.2 


0.5 


0.3 


B 


-0.5 


0.2 


0.5 


0.3 


C 


0.5 


0.2 


0.4 


0.7 


D 


1.6 


1 


0.4 


0.7 


E 


-0.7 


1 


0.4 


0.7 



TABLE I. Parameters for the different models considered. We have set j3 — 1 eV in all cases, and rj = 0.1, u>c = if the 
coupling to the dissipative bath is taken into account. Furthermore, the center of the conduction band is equal to Fermi energy 
/i L /r = ±tjV = £f for the left and right lead respectively 
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FIG. 1. Schematic representation of the tight-binding model used to parameterize the Hamiltonian, Eq. (1). The circles 
depict the atomic sites for the leads and the molecular orbital for the bridge (with energy written below) and the lines indicate 
the nonzero hopping amplitudes (written above) . 




FIG. 2. Transmission probabilities for Model A at zero bias voltage. The results shown have been obtained at different levels 
of theory: purely elastic transmittance (dotted line), numerically exact inelastic transmittance (thick dashed line), inelastic 
transmittance in the wide-band approximation (thin dashed line). The full line depicts the inelastic transmission probability 
including vibrational relaxation. 
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FIG. 3. Transmission probability for Model A as in Fig. 2. The different levels of approximation for the treatment of bath 
are shown together with exact result (solid line) and the result obtained without coupling to the bath (thick dashed line) . The 
result including up to one quantum in the final state of the bath (t^ + t^) is indistinguishable from the exact result (full t). 
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FIG. 4. Current-voltage characteristic for Model A obtained at different level of theory. 
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FIG. 5. Transmission probabilities for Model B at zero bias voltage. Shown are results obtained at different levels of theory 
as explained in the caption of Fig. 2. 
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FIG. 6. Current-voltage characteristic for Model B. Shown are results obtained at different levels of theory. 
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FIG. 7. Transmission probabilities for Model C at zero bias voltage. Shown are results obtained at different levels of theory. 
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FIG. 8. Transmission probabilities for Model C at zero bias voltage. Shown are the numerically exact results for the inelastic 
transmittance with (thick dashed line) and without (full line) vibrational relaxation, as well as the convergence of the expansion 
with respect to the number of quanta in the final state of the bath (thin lines). 
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FIG. 9. Current-voltage characteristic for Model C. Shown are results obtained at different levels of theory as explained in 
the legend. 
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FIG. 10. Current-voltage characteristic for Model C. Shown are results that illustrate the convergence of the expansion with 
respect to the number of quanta in the final state of the bath, as well as the result without vibrational relaxation (thick dashed 
line) 




FIG. 11. Transmission probabilities for the case of a bridge strongly coupled to the leads (Model D) at zero bias voltage. 
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FIG. 12. Current-voltage characteristic for Model D. 




FIG. 13. Transmission probabilities for the case of a bridge strongly coupled to the leads (Model E) at zero bias voltage. 
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FIG. 14. Potential-energy curves for the interpretation of the vibrational structures in Model E (as explained in the text). 
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FIG. 15. Current-voltage characteristic for Model E 
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